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We consider simple hydrodynamical models of galactic dark matter in which the galactic halo 
is a self-gravitating and self-interacting gas that dominates the dynamics of the galaxy. Modeling 
this halo as a sphericaly symmetric and static perfect fluid satisfying the field equations of General 
Relativity, visible barionic matter can be treated as "test particles" in the geometry of this field. 
We show that the assumption of an empirical "universal rotation curve" that fits a wide variety 
of galaxies is compatible, under suitable approximations, with state variables characteristic of a 
non-relativistic Maxwell-Boltzmann gas that becomes an isothermal sphere in the Newtonian limit. 
Consistency criteria lead to a minimal bound for particle masses in the range 30 eV < m < 60 eV 
and to a constraint between the central temperature and the particles mass. The allowed mass range 
includes popular supersymmetric particle candidates, such as the neutralino, axino and gravitino, 
as well as lighter particles (m ~ keV) proposed by numerical N-body simulations associated with 
self-interactive "cold" and "warm" dark matter structure formation theories. 
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I. INTRODUCTION 

The presence of large amounts of dark matter at the 
galactic lengthscale is already an established fact. It is 
currently thought that this dark matter is made of relique 
self-gravitating gases which are labeled as "cold" (CDM) 
or "hot" (HDM), depending on the relativistic or non- 
relativistic nature of the particles energetic spectrum at 
their decoupling from the cosmic mixture Q, [^|, ||. 
HDM scenarios are not favoured, as they seem to be in- 
compatible with current theories of structure formation 
||, ||, @- CDM, usualy examined within a newtonian 
framework, can be considered as non-interactive (colli- 
sionless particles) or self- interactive [^). N-body numeri- 
cal simmulations are often used for modeling CDM gases 
®' 0L' & &: However, in recent numerical simulations 
(see 0, |Q, |||) non- interactive CDM models present the 
following discrepancies with observations at the galactic 



scale 10 1, |llj: (a) the "substructure problem" related 
to excess clustering on sub-galactic scales, (b) the "cusp 
problem" characterized by a monotonic increase of den- 



* Contribution to Proceeedings of the IV Taller de la Division 
de Gravitacion y Fisica Matemdtica de la Sociedad Mexicana de 
Fisica (DGyFM-SMF ). Chapala, Jalisco, Mexico November 25 - 
November 30, 2001. This paper is based on the Ref. (29). 
^Electronic address: uis@nuclecu.unam.mx 

iAlso at Center for Gravitational Physics and Geometry, Penn 
State University University Park, PA, IfiSO?, ITS A.: Electronic ad- 



dress: iuncz@nuclccu.unam.mx. nuncz@gTavity.phys.p6U.edu 



§ Electronic address: gjissma i) 
^Electronic address: Itmatos' 



Q)nuclecu.uiiam.m x 
fis . cinvestav. mxj " 



sity towards the center of halos, leading to excessively 
concentrated cores. In order to deal with these prob- 
lems, the possibility of self-interactive dark matter has 
been considered, so that nonzero pressure or thermal ef- 
fects can emerge, thus leading to sclf-intcractive mod- 
els of CDM (i.e. SCDM) ||, ||, B, II, © and 
"warm" dark matter (WDM) models |17|-[g2[ that chal- 
lenges the duality CDM vs. HDM. Other proposed dark 
matter sources consist repla cing the gas of particles ap- 
proach by scalar fields j23|, [|24| and even more "exotic" 
sources p5|. 

Whether based on SCDM or WDM, current theories 
of structure formation point towards dark matter char- 
acterized by particles having a mass of the order of at 
least keV's (see Q-Q), thus suggesting that massive 
but light particles, such as electron neutrinos and axions 
(see Table 1), should be eliminated as primary dark mat- 
ter candidates (though there is no reason to assume that 
these particles would be absent in galactic halos) . Of all 
possible weakly interactive massive particles (WIMPS), 
complying with the required mass value of relique gases, 
only the massive Neutrinos (the muon or tau neutrinos), 
have been detected, whereas other WIMPS (neutralino, 
gravitino, photino, sterile neutrino, axino, etc.) are spec- 
ulative. See |§, Q H§ and Table 1 for a list of candi- 
date particles and appropriate references. 

In this paper we develop an alternative description 
of galactic DM. Since the dark matter halo constitutes 
about 90 % of the galactic mass, we consider the galac- 
tic garvitational field as a spacetime whose sole, self- 
gravitating, source is this halo, described as a perfect 
fluid. Assuming this galactic spacetime to be static and 
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sphericaly symmetric, the barionic dark matter can be 
though of as test particles following stable and circular 
geodesies of this spacetime curvature. Since the tangen- 
tial velocity, v, along these geodesies can be calculated 
for such a spacetime, we can express the field equations 
in terms of of v. However, the velocity profile v = v(r) 
has been observed for a wide range of galaxies, leading 
to "Universal Rotatition Curves" (URC's) that provide 
an empiric fit to these rotation velocities |pl|. Therefore, 
by inserting the empiric formula for the URC derived 
by Persic and Salucci Q into the field equations (given 
in terms of v(r)), we obtain a constraint on one of the 
metric coefficients. Once this constraint is solved, we 
can obtain the state variables characterizing the galactic 
dark fluid associated with this URC. We solve this con- 
straint assuming that the velocities are non-relativistic, 
thus expanding around vo/c <C 1 , where vq is the termi- 
nal velocity associated with the flattening of the URC. 



II. FIELD EQUATIONS 

Considering the line element of a static spherically 
symmetric space time 



ds 2 



-A 2 {r)c 2 dt 2 + 
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l-2M(r)/r (i) 

+r 2 (d9 2 + sin 2 9 dip 2 ) , 

the tangential velocity of test particles along stable cir- 
cular geodesic orbits can be expressed in terms of the 
metric coefficients as 



V 2 a , , rA' 
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Becomes a dynamical variable replacing A(r). Assuming 
as source of (|l|) a perfect fluid momentum energy tensor: 
T ab = (p + p)u a u b + pg ab , with u a = A- 1 S a ct , the 
following field equations in terms of (0) become 
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where k = 8ttG /c 4 and a prime denotes derivative with 
respect to r. Writing the field equations in terms of the 
orbital velocity, v, provides a useful insight into how an 
(in principle) observable quantity relates to spacetime 
curvature and with physical quantities (state variables) 
which characterize the source of spacetime. 



III. THERMODYNAMICS 

If we assume that the self gravitating ideal "dark" gas 
exists in physical conditions far from those in which the 
quantum properties of the gas particles are relevant, we 
would be demanding that these particles comply with 
Maxwell-Boltzmann (MB) statistics. Following p4| , the 
condition that justifies an MB distribution is given by 



n h 



(mk B T) 3 / 2 



< 1 
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where n, T, h and k B are, respectively, the particle num- 
ber density, absolute temperature, Planck's and Boltz- 
mann's constants. If the constraint @) holds and we 
further assume thermodynamical equilibrium and non- 
relativistic conditions, the ideal dark gas must satisfy the 
equation of state of a non-relativistic monatomic ideal gas 



p = mc n + —nk B T 



P : 



(7) 



whose macroscopic state variables can be obtained from 
a MB distribution function under an equilibrium Kinetic 
theory approach (the non-relativistic and non-degenerate 
limit of the Jiittner distribution) [ p2[ . An equilibrium 
MB distribution restricts the geometry of spacetime |?3| , 
resulting in the existence of a timclike Killing vector field 
(3 a = (3u a , where /3 = mc 2 /k B T, as well as the following 
relation (Tolman's law) between the 4-acceleration and 
the temperature gradient 



ii a + h b a (lnT), b = 
leading to 
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The particle number density n trivialy satisfies the con- 
servation law J a . a = where J a = n u a , thus the number 
of dark particles is conserved. Notice that given (|J) and 
(||) , the equation of state (0) and the temperature from 
the Tolman law (ph , we have two different expressions for 
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-8 vv' r + (v 2 + 9) (2 v 2 + l)] M 
K mc 2 r 3 (1 + v 2 ) 

+4vv' r 2 - v 2 (7 + v 2 ) r 
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The quantity mc 2 n in ( |ll| ) follows directly from equa 
tions (| ) and (||) , while n in (|To|) also follows from p in 



with A oc exp[J (w 2 /r)dr]. Consistency requires that 
ro|) and (O) yield the same expression for n. 



IV. DARK FLUID HYDRODYNAMICS 



We shall assume for v 2 the empiric dark halo rotation 
velocity law given by Persic and Salucci , |3l[] 



,,2 J2. 
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where r opt is the "optical radius" containing 83 % of 
the galactic luminosity, whereas the empiric parameters 
a and vq, respectively, the ratio of "halo core radius" 
to r opt and the "terminal" rotation velocity, depend on 
the galactic luminosity. For spiral galaxies we have: 
Vq = v 2 pt (l — + a 2 ) , where v opt — v(r opt ) and the best 
fit to rotation curves is obtained for: a = 1.5 (L/L*) 1 / 5 
and (3 = 0.72 + 0.44 log 10 (£/£*), wh ere L* = 1O 1O - 4 L . 
The range of these parameters for spiral galaxies is 
125 km/sec < vq < 250 km/sec and 0.6 < a < 2.3. 
Inserting @ into (|) and (|) we obtain 
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FIG. 1: Comparison of n obtained from \1<\) and l{2U\). This 
plot displays the adimensional quantity Kmc r„ p n, as a func- 
tion of x, obtained from truncating the right hand sides of 
JT^ j (solid curve) and (dotted curve with crosses) up to 
fourth order in vo/c and assuming the consistency condition 
Since we are using an empiric law for observed galactic 
rotation curves (the URC given by ^1^)), the fact that these 
two expressions for n are so close to each other provides an 
empirical justification for the compatibility between MB dis- 
tribution and these rotation curves. 
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we obtain the expanded forms of p and p by inserting (|lz 
and ( |l4| ) into (^) and (||) and then expanding around v , 
leading to 
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where T c = T(0) and we have set an integration constant 
to zero in order to comply with the consistency require- 
ment that vq = implies flat spacetime (A = 1, M = 0). 
Since the velocities of rotation curves are newtonian, 
!)o < c (typical values are v /c « 0.5 x 10~ 3 ), instead 
of evaluating dl4| ) we will expand this quadrature around 
vq/c (in order to keep the notation simple, we write t>o 
instead of vq/c). This yields 
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while the expanded form for T follows from (Fti 
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FIG. 2: Minimal mass for which the Maxwell-Boltzmann dis- 
tribution is applicable. This graph displays m (in keV's) 
as a function of vq/c and b = ax, respectively, the termi- 
nal velocity and 'halo core radius' associated with the URC 
given in @. Assuming typical ranges for sptral galaxies: 
125 km/sec < vo < 250 km/sec and 1 kpc < b < 5 kpc, we ob- 
tain masses in the range of 30 eV < m < 60 eV that follow 
from the right hand side of the relation \2^), providing the cri- 
terion for applicability of the Maxwell-Bolzmann distribution. 
Dark matter particle candidates complying with an MB dis- 
tribution must have much larger mass than the plotted values 
30eV<m< 60 eV. 
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In order to compare n obtained from ( |10[ ) and (Jllj), we 
substitute ( p"2| ) and ([u]) into (|ll|) and expand around Vq, 
leading to 
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Since vq/c -C 1, a reasonable approximation is ob- 
tained if the leading terms of n from ([l9]) and (20) co- 
incide. By looking at these equations, it is evident that 
this consistency requirement implies 
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where vq denotes a velocity (cm/sec) and not the adi- 
mensional ratio Vq/c. Since higher order terms in vq/c 
have a minor contribution, the two forms of n are ap- 
proximately equal. This is shown in Figure 1 displa ying 
the adimensional quantity k mc 2 n r 2 pt from ( |l9| ) and (|20f ) 
as functions of x for typical values Vq/c = 0.0006, a = 1 
and eliminating T c with (|2l]). Equation ( |l4| ) shows how 
"flattened" rotation curves, as obtained from the empiric 
form (|lj), lead to M oc r 3 for r « and M oc r for large 
r. Equations @ to @ represent a relativistic gener- 
alization of the "isothermal sphere" that follows as the 
newtonian limit of an ideal Maxwell-Boltzmann gas char- 
acterized by p s» mc 2 n,p <C p and T w T c . In fact, using 
newtonian hydrodynamics we would have obtained only 
the leading terms of equations ( |l5| ) to (pil|). It is still 
interesting to find out that the isothermal sphere can be 
obtained from General Relativity in the limit Vq/c <C 1 
by demanding that rotation curves have a form like jl^). 
The total mass of the galactic halo, usualy given as M 
evaluated at the radius r = r 2 oo (the radius at which p is 
200 times the mean cosmic density) . Assuming this den- 
sity to be Ri 10~ 29 gm/cm 3 together with typical values 
vq = 200km/sec and a = 1 yields r2oo ~ 150 kpc. Eval- 
uating M at this values yields about 10 17 M Q , while M 
evaluated at a typical "optical radius" r — 15 kpc leads 
to about 10 12 Mq, an order of magnitude larger than the 
galactic mass due to visible matter. 
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V. DISCUSSION. 

So far we have found a reasonable approximation for 

(19) 

galactic dark matter to be described by a self gravitat- 
ing Maxwell-Boltzmann gas, under the assumption of the 
empiric rotation velocity law (|l|). The following consis- 
tency relations emerge from equations (jll]) , fl2C| ) and ( pl| ) 



while n in ( pX) f ollows by substituting ( |l4| ) into (17), 
using T from (13) and then expanding around vq. This 
yields 
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TABLE I: Particle candidates for a MB Dark Matter gas. 



SCDM/WDM 
Light Candidates 
Light Gravitino 

Sterile Neutrino 



Standard Neutrinos 
Light Dilaton 
Light Axino 
Majoron 

Mirror Neutrinos 



CDM 
Heavy Candidates 
Neutralino 

Axino 
Gravitino 




hence, bearing in mind that n < n c and T w T c , the con- 



dition ( 
with (t 



3) for the validity of the MB distribution together 
I) yields the condition 



m » 
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a criteria of aplicability of the MB distribution that is 
entirely given in terms of m, the fundamental constants 
G, h and the empiric parameters vq and a r opt (the "ter- 
minal" rotation velocity and the "core radius" ) [^9| . For 
dark matter dominated galaxies (spiral and low surface 
brightness (LSB)) [Q these parameters have a small 
variation range: r opt « 15kpc, 0.6 < a < 2.3 and 
125 km/sec < vq < 300 km/sec, the constraint (|23|) does 
provide a tight estimate of the minimal value for the mass 
of the particles under the assumption that these particles 
form a self gravitating ideal dark gas complying with MB 
statistics. As shown in Figure 2, this minimal value lies 
between 30 and 60 eV, thus implying that appropriate 
particle candidates must have a much larger mass than 
this figure [ p9| . This minimal bound excludes, for in- 
stance, light mass thermal particles such as the electron 
neutrino (m Ve < 2.2 eV). The axion is also very light 
{m,A ~ 10~ 5 eV) but it is not a thermal relique and 
so we cannot study it under the present framework. The 
currently accepted estimations of cosmological bounds on 
the sum of masses for the three active neutrino species is 
about 24 eV, a value that would apparently rule out all 
neutrino flavours. However, recent estimations of these 
cosmological bounds have raised this sum to about 1 keV 
|j42f , hence more massive neutrinos could also be accomo- 
dated as dark matter particle candidates. Estimates of 
masses of various particle candidates are displayed in Ta- 
ble 1. 



FIG. 3: Relation between particle mass and central temper- 
ature. This graph displays the relation between log ln (T c ) fin 
K) and log 10 (m) (in keV's) that follows from equation fyzjj ) 
for a terminal velocity vq = 200 km/sec. Almost identical 
plots are obtained for other velocities in the observed range 
125 km/sec < vo < 250 km/sec. The circle and box symbols re- 
spectively denote the proton and electron mass yielding central 
temperatures of the order T c « 10 6 , 10 3 K. The central tem- 
perature for light particles in the range 0.5 keV < m < 100 keV 
is less than 300 K (rectangle in the left), while for massive 
sypersymmetric particles in the range 1 GeV < m < 100 GeV, 
we have T c as large as 10 9 K (rectangle on the right). How- 
ever, such high temperatures cannot rule out these weakly in- 
teractive particles as components of the dark matter MB gas. 



Since T w T c , the consistency condition ( pl| ) provides 
the following constraint on the temperature and particles 
mass of the dark gas 
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(24) 



where we have taken vq = 300 km/sec. Considering in 
( p4| ) the minimal mass range that follows from (|2^), we 
would obtain gas temperatures consistent with the as- 
sumed typical temperatures of relic gases: T c w 2 to 4K. 
However, since we have no way of infering a value for the 
temperature of the ideal dark gas, we have no clear cut 
criterion for the estimate of a maximal bound for this 
mass. If we assume that the ideal dark gas is made of 
electrons or barions, so that m — m r , or m = m e , then 



condition (23) for applicability of the MB distribution is 
certainly satisfied and ( pMj ) implies a temperature of the 
order of T c « 10 3 K for electrons and T c w 10 6 K for 
barions. Obviuosly, barions or electrons at such a high 
temperatures would radiate and certainly not remain un- 
observably "dark" . However, as long as the interaction 
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is weak and the particles are not charged, we cannot rule 
out any other particle candidate only on the basis of the 
gas temperature, even if this temperature is very high 
(see Figure 3) as in the case of massive supersymmetric 
particles. As shown in Table 1, a wide range of weakly 
interactive particles can be considered as possible main 
components of a MB dark gas, including popular super- 
symmetric particles (the neutralino), as well as hypothet- 
ical light particles predicted by current literature based 
on WDM models of structure formation Q , ||| . The main 
novelty of the present paper is the fact that it is based on 
a general relativistic hydrodynamics, as opposed to nu- 
merical simulations , newtonian or Kinetic Theory 
perturbative approaches (see [Ti^-p^]). 

Finaly, the fact that we have obtained a minimal mass 
on the range 30 — 60 eV , that seems to discriminate 
against very light thermal particles like the electron neu- 
trino, coincides with the fact that these HDM particle 
candidates tend to be ruled out because of their inabil- 
ity to produce sufficient matter clustering 0, ||, p9| . 
In spite of these arguments, a self gravitating gas of this 



type of particles accounting for a galactic halo, would 
have to be modeled, either as a relativistic MB gas (very 
light particles can be relativistic even at low tempera- 
tures) and/or in terms of a distribution that takes into 
account Fermi-Dirac or Bose-Einstein statistics. These 
studies will be undertaken in future papers [^9| . 
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